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A fully relativistic finite nucleus wave-function approach to the non-mesonic weak decay of the 
a 2 C hypernucleus is presented. The model is based on the calculation of the amplitudes of the 
tree-level Feynman diagrams for the AN — > NN process and includes one-pion exchange and one- 
kaon exchange diagrams. The pseudo-scalar and pseudo-vector choices for the vertex structure 
are compared. Final-state interactions between each one of the outgoing nucleons and the residual 
nucleus are accounted for by a complex phenomenological optical potential. Initial AiV and final NN 
short-range correlations are included by means of phenomenological correlation functions. Numerical 
results are presented and discussed for the total non-mesonic decay width T nm = r n +T p , the T n /T p 
ratio, the oa intrinsic asymmetry parameter, and the kinetic energy and angular spectra. 

PACS numbers: 21.80.+a Hypernuclei; 24.10. Jv Relativistic Models 

I. INTRODUCTION 

The birth of hypernuclear physics dates back to 1952 [1] when the first hypernuclear fragment originated from 
the collision of a high-energy cosmic proton and a nucleus of the photographic emulsion exposed to cosmic rays was 
observed through its weak decays, revealing the presence of an unstable particle: this was interpreted as due to 
the formation of a nucleus in which a neutron is replaced by the A hyperon, i.e., the lightest strange baryon. A 
hypernucleus is a bound system of neutrons, protons, and one or more hyperons. Only the lightest hyperon, the A, 
is stable with respect to esoenergetic strong and electromagnetic processes in nuclear systems. Therefore, the most 
stable hypernuclei are those made up of nucleons and a A particle. We denote with a hypernucleus with Z 
protons, (A — Z — 1) neutrons, and a A (A- hypernucleus). 

Hypernuclei represent a unique laboratory for the study of strong and weak interactions of hyperons and nucleons 
through the investigation of hypernuclear structure and decay. The A particle, which does not have to obey the 
Pauli principle, is an ideal low-energy probe of the nuclear environment which allows a deepening of classical nuclear 
physics subjects, such as the role of nuclear shell models and the dynamical origin of the nuclear spin-orbit interaction. 
Hypernuclear physics also establishes a bridge between nuclear and hadronic physics, since many related issues can in 
principle unravel the role played by quarks and gluons partonic degrees of freedom inside nuclei. In this direction, the 
study of hybrid theories combining meson-exchange mechanisms with direct quark interactions have the potentiality 
to teach us something on the confinement phenomenon, an issue still far from being satisfactorily understood. 

In A-hypernuclei the A can decay via either a mesonic or a non-mesonic strangeness-changing weak interaction 
process. In the nuclear medium the mesonic decay, A — > Ntt, which is the same decay of a free A, is strongly 
suppressed, but in the lightest hypernuclei, by the effect of the Pauli principle on the produced nucleon, whose 
momentum (~ 100 MeV/c) is well below the Fermi momentum. In the non-mesonic weak decay (NMWD) the pion 
produced in the weak A — » Ntt transition is virtual and gets absorbed by neighbor nucleons. Then, two or three 
nucleons with high momenta (~ 400 MeV/c) are emitted. We can distinguish between one and two-nucleon induced 
decays, according to whether the A interacts with a single nucleon, either a proton, Ap — > np (decay width r p ), or a 
neutron, An — > nn (T n ), or with a pair of correlated nucleons, ANN — > nNN (T2). Mesons heavier than the pion 
can also mediate these transitions. The NMWD process is only possible in the nuclear environment and represents 
the dominant decay channel in hypernuclei beyond the s-shell. The total weak decay rate is given by the sum of the 
mesonic (r m ) and non-mesonic (r nm ) contributions: 

^tot — r. m -f r nm , (i) 

with 

r„ m = r 1 + r 2 , r 1 = r p + r„. (2) 

The fundamental interest in the NMWD mode is that it provides a unique tool to study the weak strangeness 
changing (\A.S ~ 1|) baryon-baryon interaction AA — > nN, in particular its parity conserving part, that is much more 
difficult to study with the weak JVJV — > NN transition, that is overwhelmed by the parity-conserving strong JVJV 
interaction. Since no stable hyperon beams are available, the weak process AA — > nN can be investigated only with 
bound strange systems. The study of the inverse process pn — > pA would however be useful. 
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Although the relevance of the NMWD channel was recognized since the early days of hypernuclear physics, only in 
recent years the field has experienced great advances due to the conception and realization of innovative experiments 
and to the development of elaborated theoretical models [2-8]. 

For many years the main open problem in the decay of hypernuclei has been the T n /T p puzzle, i.e., the disagreement 
between theoretical predictions and experimental results of the ratio between the neutron- and proton-induced decay 
widths: for all the considered hypernuclei the experimental ratio, in the range ~ 0.5-^2, was strongly underestimated 
(by about one order of magnitude) by the theoretical results. The T n /T p ratio directly depends on the isospin 
structure of the weak process driving the hypernuclear decay. The analysis of the ratio is a complicated task, due 
to difficulties in the experimental extractions, which require the detection of the decay products, especially neutrons, 
and to the presence of additional competing effects, such as final-state interactions (FSI) of the outgoing nuclcons 
and two-nucleon induced decays, which could in part mask and modify the original information. 

In the first theoretical calculations the one-pion-exchange (OPE) nonrelativistic picture was adopted as a natural 
starting point in the description of the AA — > nN process, mainly on the basis of its success in predicting the basic 
features of the strong NN interaction. The first OPE models were able to reproduce the non-mesonic decay width 
r„ m = r„ + T p but predicted too small r„/r p ratios [9-14]. It thus seemed that the theoretical approaches tend 
to underestimate T n and overestimate T p . A solution of the puzzle then requires devising dynamical effects able to 
increase the n-induced channel and decrease the p-induced one. 

In the following years the theoretical framework was improved including the exchange of all the pseudo-scalar and 
vector mesons, in the form of a full one-meson-exchange (OME) model, or properly simulating additional effects, above 
all initial short-range correlations (SRC) and FSI, by means of direct quark mechanisms and many-body techniques 
[9, 10, 12, 13, 15-24]. In particular, the inclusion of K exchanges seems essential to improve the agreement between 
theory and experiments. Only a few of these calculations have been able to predict a sizeable increase of the T n / T p 
ratio [13-15, 24], but no fundamental progress has been achieved concerning the deep dynamical origin of the puzzle. 

The situation has been considerably clarified during the very last years, thanks to considerable progress in both 
experimental techniques [25-31] and theoretical treatments [13-15, 24, 32-40]. From the experimental point of view, 
the new generation of KEK experiments has been able to measure the fundamental observables for the \He and 
hypernuclei with much more precision as compared with the "old" data, also providing the first results of simultaneous 
one-proton and one-neutron energy spectra, which can be directly compared with model calculations. Very recently, 
it has also been possible to obtain for the first time coincidence measurements of the nucleon pairs emitted in the 
non-mesonic decay, with valuable information on the corresponding angular and energy correlations. These new data 
further refine our experimental knowledge of the hypernuclear decay rates, also allowing a cleaner and more reliable 
extraction of the r„/r p ratio. From the theoretical point of view, crucial steps towards the solution of the puzzle 
have been carried out, mainly through a non-trivial reanalysis of the pure experimental results by means of a proper 
consideration of FSI and rescattering mechanisms, inside the nuclear medium, for the outgoing nucleons, as well 
as of the two-nucleon induced channel. This strict interplay between theory and experiments is at the basis of the 
present belief that the Y n / Y p puzzle has been solved. In particular, this is due to the study of nucleon coincidence 
observables, recently measured at KEK [29, 30], whose weak-decay-model independent analysis carried out in [32, 33] 
yields values of T n / T p around 0.3 0.4 for the \He and \ 2 C hypernuclei, in satisfactory agreement with the most 
recent theoretical evaluations [13-15, 24]. New, more precise results are expected from forthcoming experiments at 
DA$NE [41] and J-PARC [42]. 

Another intriguing issue is represented by the asymmetry of the angular emission of non-mesonic decay protons from 
polarized hypernuclei. The large momentum transfer involved in the n(jr + , A + )A reaction can be exploited to produce 
final hypernuclear states characterized by a relevant amount of spin-polarization, preferentially aligned along the axis 
normal to the reaction plane [43, 44]. The hypernuclear polarization mainly descends from a non- negligible spin- flip 
term in the elementary mr + — > AK + scattering process, which in turn interferes with the spin-non-flip contribution 
[45]: in free space, and for \p n \ = 1.05 GcV and 9k — 15°, the final hyperon polarization is about 75%. 

Polarization observables represent a natural playground to test the present knowledge of the NMWD reaction 
mechanism, being strictly related to the spin-parity structure of the elementary AA — > nN interaction. Indeed, by 
focusing on the p-induced channel, experiments with polarized hypernuclei revealed the existence of an asymmetry 
in the angular distributions of the emitted protons with respect to the hypernuclear polarization direction. Such an 
asymmetry originates from an interference effect between parity-violating and parity-conserving amplitudes for the 
Ap — > np elementary process, and can thus complement the experimental information on the T n and T p partial decay 
rates, which are instead mainly determined by the parity-conserving contributions. As for the T n / T p ratio, FSI could 
play a crucial role in determining the measured value of this observable. 

The asymmetry puzzle concerns the strong disagreement between theoretical predictions and experimental extrac- 
tions of the so-called intrinsic asymmetry parameter a\. The first asymmetry measurements [43, 44] with limited 
statistics gave large uncertainties and even inconsistent results. The very recent and more accurate data from KEK- 
E508 [5, 47, 48] favour small values of a\, compatible with a vanishing value. Moreover, the observed asymmetry 
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parameters are negative for A 2 C and positive (and smaller, in absolute value) for \He. Theoretical models generally 
predict negative and larger values of a\. FSI effects do not improve the agreement with data [49]. The inclusion, 
within the usual framework of nonrelativistic OME models, of the exchange of correlated and uncorrelated pion pairs 
[50] greatly improves the situation. Indeed, it only slightly modifies the non-mesonic decay rates and the T n / T p ratio, 
but the modification in the strength and sign of some relevant decay amplitudes is crucial and yields asymmetry 
parameters which lie well within the experimental observations. In particular, a small and positive value is now 
predicted for a\ in A He. This important achievement justifies the claim that also the asymmetry puzzle has finally 
found a solution. 

Recent experimental and theoretical studies have led to a deeper understanding of some fundamental aspects of 
the NMWD of A-hypernuclei. From a theoretical point of view, the standard approach towards these topics has been 
strictly nonrelativistic, with both nuclear matter and finite nuclei calculations converging towards similar conclusions: 
nonrelativistic full one-mcson-exchange plus two-pion-exchange models, based on the polarization-propagator method 
(PPM) [51, 52] or on the wave-function method (WFM) [11, 13, 14], seem able to reproduce all the relevant observables 
for the \He and A 2 C light- medium hypernuclei. A crucial contribution to this achievement is however due to a non- 
trivial theoretical analysis of KEK most recent coincidence data, based on the proper consideration and simulation 
of nuclear FSI and two-nucleon induced decays [33, 50, 53, 54]. In those nonrelativistic models many theoretical 
ingredients are included with unavoidable approximations. Initial-state interactions and strong AN SRC are treated 
in a phenomenological way, though based on microscopic models calculations. The inclusion of the full pseudo-scalar 
and vector mesons spectra, in particular the strange K and K* mesons, in the context of complete OME nonrelativistic 
calculations, somewhat clashes with the still poor knowledge of the weak A- TV-meson coupling constants for mesons 
heavier than the pion. Their evaluation requires model calculations which unavoidably introduce a certain degree of 
uncertainty in the corresponding conclusions about the r„/r p ratio and other observables. In order to reproduce few 
observables, i.e., the decay rates T nm , T n , and T p , and the asymmetry parameter a A for \He and \ 2 C, these models 
need to include many dynamical effects, such as the exchange of all the possible mesons, plus two-pion exchange, plus 
phenomenological a mesons, plus the corresponding interferences, and, moreover, strong nuclear medium effects in 
the form of non-trivial FSI. Although the inclusion of many theoretical ingredients can be considered as a natural and 
desirable refinement of the simple OPE models, all the improvements do not seem to provide a significantly deeper 
insight into the decay dynamics. Despite all the theoretical efforts, the solution of the r„/r p and cla parameter 
puzzles seems due to effects, such as distortion, scattering and absorption of the primary nucleons by the surrounding 
nuclear medium, rather than to the weak-strong interactions driving the elementary AN — > nN or ANN — > nNN 
process. The re-analysis of the recent KEK experimental data [25, 29, 30, 55] and the corresponding extraction of 
r n /r p are indeed completely independent of the weak- decay- mechanisms [33, 38], but depend strongly on the model 
adopted to describe FSI and on somewhat arbitrary assumptions, e.g. on the ratio r 2 /Fi between two-nucleon and 
one-nucleon induced non-mesonic decay rates. 

The presently available experimental information on hypcrnuclcar decay is still limited and affected by uncertainties 
of both experimental and theoretical nature. Moreover, single-nucleon spectra seem to point at a possible systematic 
protons underestimation [38] . The new generation of experiments planned in various laboratories worldwide is expected 
to produce more precise data on the already studied observables as well as new valuable information in the form of 
differential energy and angular decay particles spectra. 

In spite of the recent important achievements, the NMWD of hypernuclei deserves further experimental and the- 
oretical investigation. From the theoretical point of view, the role of relativity is almost unexplored. But for a few 
calculations in [46, 56, 57, 64] no fully relativistic model has been exploited to draw definite conclusions about the 
role of relativity in the description of the weak decay dynamics. 

In this paper we present a fully relativistic model for the NMWD of \ 2 C [58]. The adopted framework consists of 
a finite nucleus WFM approach based on Dirac phenomenology. As a first step the model includes only OPE and 
one-kaon-exchange (OKE) diagrams, and is limited to one-nucleon induced decay. We are aware that the neglected 
contributions could play an important role in the decays. Our aim is to explain all the at least qualitative features 
of the hypernuclear NMWD with a conceptually simple model, in terms of a few physical mechanisms and free 
parameters. We stress that, dealing with a fully relativistic treatment of the weak dynamics based on the calculation 
of Feynman diagrams within a covariant formalism, it is quite difficult to directly compare such an approach and its 
results to standard nonrelativistic OME calculations. We will thus rather focus on the internal coherence and on the 
theoretical motivations of the model. 

The model is presented in Sec. 2. Numerical results for the total non-mesonic decay width T nm , the r„/r p ratio, 
the a \ intrinsic asymmetry parameter, as well as for kinetic energy and angular spectra are presented and discussed 
in Sec. 3. The sensitivity to the choice of the main theoretical ingredients is investigated. The theoretical predictions 
of the model are compared with the most recent experimental results. Some conclusions are drawn in Sec. 4. 
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II. MODEL 

In this Section we present a fully relativistic finite nucleus wave-function approach to study the NMWD of the 
hypernucleus. Our model is based on a fully relativistic evaluation of the elementary amplitude for the A./V — > NN 
process, which, at least in the impulse approximation, is the fundamental interaction responsible for the NMWD. 
Covariant, complex amplitudes are calculated in terms of proper Feynman diagrams. The tree-level diagram involves 
a weak and a strong current, connected by the exchange of a single virtual meson. Integrations over the spatial 
positions of the two vertices as well as over the transferred 3-momentum are performed. As a first approximation, 
only OPE and OKE diagrams are considered. Possible two-nucleon induced contributions are neglected, even if they 
could play an important role in the hypernuclear decay phenomenology. 

Interested readers can find further details about the present model in the PhD thesis of Ref. [58] , where an extensive 
analysis of the adopted formalism as well as of the involved theoretical ingredients is provided. 

In the calculation of the hypernuclear decay rate the Feynman amplitude must then be properly included into a 
many-body treatment for nuclear structure. The amplitude is therefore only a part of the complete calculation, but 
it is the basic ingredient of the model and involves all the relevant information on the dynamical mechanisms driving 
the decay process. 

Short range correlations are also included in the model, coherently with what commonly done in most nonrelativistic 
calculations, since the relatively high nucleon energies involved in the hypernuclear NMWD can in principle probe 
quite small baryon-baryon distances, where strong interactions may be active and play an important role. Following 
a phenomenological approach, we have chosen to include initial SRC effects by means of a multiplicative local and 
energy-independent function, whose general form [46] provides an excellent parametrization of a realistic AiV cor- 
relation function obtained from a G-matrix nonrelativistic calculation [59, 60]. The problem of ensuring a correct 
implementation of such a nonrelativistic SRC function within a relativistic, covariant formalism has been addressed 
in Ref. [61] and shown to be tightly connected with the choice of the interaction vertices. For full generality, we also 
choose to account for possible strong NN short range interactions acting on the two final emitted nucleons, again 
adopting a simple phenomenological average correlation function [13] which provides a good description of nucleon 
pairs in 4 Hc [62] as calculated with the Reid soft-core interaction [63]; such final-state SRC could in principle play 
an important role, and they complement the final-state interactions between each of one of the two emitted nucleons 
and the residual nulceus, that is accounted for in our model by a relativistic complex optical potential. 

A. Coupling ambiguities 

In order to devise a relativistic treatment of the elementary AA" — > NN process, great care must be devoted 
to the choice of the Dirac-Lorentz structure for the strong and weak parity-conserving vertices. The pseudo-scalar 
(PS) prescription, that consists in a 175 Dirac structure, and the pseudo-vector (PV) one, that contains a 757^9^ 
axial-vector structure, are in principle equivalent, at least for positive energy on-shell states, because they descend 
from equivalent Lagrangians. However, ambiguities arise when one tries to take into account SRC in terms of a 
multiplicative local and energy- independent function f(r). Such ambiguities are not of dynamical origin and should 
not be mistaken as relativistic effects: they are simply bound to the phenomenological way of including (initial 
and final) short range correlations, by matching a nonrelativistic correlation function within a relativistic Feynman 
diagram approach. The crucial observation, in this regard, is that it is possible to give theoretical reasons [61] to 
prefer the PV coupling in its modified version where the 4-derivative operates on the propagator (PV'), over the 
PS coupling and also over the standard PV one, where the 4-derivative acts on the matrix element. On the one 
hand, the PV choice permits to recover, in the nonrelativistic limit, the standard OPE potential, multiplied by f(r), 
which is commonly used as the starting point in nonrelativistic calculations, whereas the PS and PV couplings yield 
a simple Yukawa function in the same limit: this allows, at least in principle, a comparison between relativistic and 
nonrelativistic results. On the other hand, a microscopic model of (initial) SRC effects, adopting standard PS NNtt 
vertices and introducing an additional w-exchange mechanism simultaneous to the OPE dominant one, produces a 
result analogous to what can be derived in a phenomenological tree-level approach contemplating the inclusion of 
a SRC function, provided in this case the modified derivative PV coupling, rather than the PS one, is employed. 
The main feature is the development of an explicit dependence of the interaction matrix elements on the exchanged 
three-momentum (through the momentum involved in the corresponding loop integrals, in the microscopic model, or 
the derivative effect of the PV coupling, in the tree-level phenomenological approach). When dealing with a more 
complex model for nuclear structure, we do not generally use positive energy on-shell states. Still the general message 
keeps its validity, though the details of the explicit calculations may be different. In order to correctly treat SRC, 
nuclear currents showing a dependence on the 3-momentum transfer q are needed, which in the simple model above 
correspond to matrix elements between external spinors and intermediate spinors carrying the momentum of the 
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intermediate state excited by the heavy meson. This can be achieved using the PV' coupling acting on the pion field, 
while the use of the PS or standard PV couplings, as done for instance in Ref. [46], would generate nuclear currents 
independent of q, corresponding, within the considered simple SRC model, to matrix elements between spinors all 
carrying the external momenta. 



B. Pseudo-scalar couplings 

As a first example, we employ a PS coupling for the strong vertex and for the parity-conserving part of the weak 
interaction. The AN — > nN fundamental process can then be decomposed into a weak ANtt vertex, governed by the 
weak Hamiltonian 

= iG F ml tt« (A n + B-k 75 ) t ■ fa , (3) 
and a strong NNtt vertex, driven by the Hamiltonian 



ft&k = 1 9NN, 75 r • fa * W . (4) 



The Dirac spinors and <J>^ ) are the wave functions of the bound A hyperon and nucleon inside the hypernucleus, 

&n is the Dirac spinor representing the scattering wave function of each one of the two final nucleons, r is the vector 
formed by the three Pauli matrices, and fa is the isovector pion field. The Fermi weak constant Gf and the pion 
mass TOtt give Gp — 2.21 x 1CT 7 . The empirical constants = 1.05 and = —7.15 are adjusted to the free A 
decay and determine the strengths of the parity-violating and parity-conserving non-mesonic weak rates, respectively. 
Finally, gNN-n = 13.16 is the strong NNtt coupling. The initial A and final nucleon fields, and , are defined 
in space-spin-isospin space and they are described by a space-spin part times a two component isospinor. In addition, 
the A field is represented as a pure m tA — — 1/2 state to enforce the il AI = 1/2 selection rule. 

The relativistic Feynman amplitude for the two-body matrix element describing the AN — > NN transition, driven 
by the exchange of a virtual pion, can be written as 

T}ff = iG F mlg NNv f d 4 x J d 4 y f™(\x - y\) 
6° b A„(x-y)fgZ(\x- y \) 

where ^ a _ 1 , 2 and ^aN,/iN,m tN are the bound A and nucleon wave functions, with ajv,A = {nlj}N,A quantum 
numbers and total spin (isospin) projections hn,a ( r nt N , mt A = —1/2), and , with i = 1,2, are the 

scattering wave functions for the two final nucleons emitted in the hypernuclear NMWD, with asymptotic momenta 
hi and spin (isospin) projections m Si (mt i ) . In both the initial and the final baryon wave functions it is possible to 
factor out the isospin 2-spinors as well as the energy-dependent exponentials: ^(a;) = ip(x) e ~ lEx xj™ ; t 1 y 2 , where E 
is the total energy of the considered baryon. The factor (\x — y\) represents a short-range two-body correlation 
function acting on the initial A and N baryons, and similarly f^(\x — y\) describes possible short range interactions 
between the two final nucleons emerging from the interaction vertex. (x — y) is the Fourier transform of the product 
of the pion propagator with the vertex form factors (supposed to be equal for the strong and weak vertices), i.e., 

After performing time integrations in Eq. (5) and taking advantage of the q° part of the integral in Eq. (6), we get 
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for the relativistic amplitude the expression 

Tjff = i G F ml g NN7I 1 J d?x J d 3 y /™(|x - y|) 

x ^{\x- v \) f f N in N {\x- v \) [^l m Jy)j 5 ^i,, N (y) 

x (2tt) <5 (£1 +E 2 -E A - E N ) , 



where 



A n (\x-y\) 



iq(x-y) 



(2-7r) 3 (g ) 2 - q 2 - ml + iz 



(7) 



(8) 



with q° — E\ — Ei = E 2 — Em and I is an isospin factor that depends on the considered decay channel (either p- or 
n- induced), i.e., 



I = 



mt N 
Xl/2 



(9) 



It is easy to check that I is different from zero only for the charge-conserving processes Ap — > np and An — > nn. For 
the calculation of the integral over the 3-momentum transfer q in Eq. (7) we choose a monopolar form factor, i.e., 



2 A^-r^ 

^ )- A2-Q2 



(10) 



where m v ~ 140 MeV is the pion mass and A ff ~ 1.3 GeV is the cut-off parameter [46]. The initial AAT correlation 
function adopted in our calculations is [46] 



m(r) = (l - e- r2 / a2 y + br 2 e- r ^ c2 , 
with n = 2, a = 0.5, b = 0.25, c = 1.28, while the final AW correlation function is chosen as [13] 

/j&(r) = l-j («cr), 



(11) 



(12) 



where jo(x) = is the first spherical Bcssel function, and q c — 3.93 fm _1 . The two correlation functions of Eqs. 
(11) and (12) are plotted in Fig. 1. As a consequance of the approximations adopted in the present calculation, we 
could then, for practical purposes, treat initial and final SRC as a whole, in terms of an overall correlation functions 
defined as 



f(\x-y\) = r A n N (\x-y\)f N i N (\x-y\). 



(13) 



T (PV') 



C. Pseudo- vector couplings 

When we use derivative PV couplings, the equivalent of Eq. (5) is 

= iG F m\ g NN7T J ' d 4 x J d 4 y 

fei,m si ,m tl v ; V 2M J a A ,MA,-l/2 v - > 

x 5 ab A n (x-y)f(\x-y\) 



(y) 



2M N 



(14) 
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FIG. 1: Initial AiV f A n ^(r) (Eq. (11), solid line) and final AW /#£(r) (Eq. (12), dashed line) correlation functions, with r 
representing the baryon-baryon relative distance. 



where M = (M N + M A )/2, with M A = 1.1156 GeV, and now f(\x - y\) is defined as in Eq. (13). In Eq. (14) the 
space-time derivatives act just on the pion propagator A 7T (x — y) (given in Eq. (8)) and not on the SRC function 
f(\x — y\), which is considered as a phenomenological ingredient entering Eq. (14) in a factorized form. Thus, the 
involved derivatives translate into multiplicative terms under the integral over q. After the time integrations in Eq. 
(14), we obtain the final expression of the relativistic amplitude 



T}?% ) = iG F m\ g NN7I 1 I d 3 x I d 3 y 



J (27T 



d 3 q e 



iq(x-y) 



) 3 q 2 - ml + ie 



r&tf) f(\x-y\ 



£>7T 



(s) 
k 2 ,r, 



(y) 



2M, 



N 



x (2tt) S (E 1 +E 2 -E a - E n ) 



(15) 



where I, ^(q 2 ), and f(\x — y\) are defined in Eqs. (9), (10), and (13), respectively. Eq. (15) must be evaluated at 
q° = q°=E A - E 1 = E 2 - E N . 

The crucial difference between Eq. (7) and Eq. (15) is that in Eq. (15), obtained adopting derivative PV couplings 
at the vertices, the matrix elements between the initial bound and the final scattering states explicitly depend on the 
3-momcntum transfer q, while using PS couplings the matrix elements in Eq. (7) are independent of q. 

Though representing a computational complication, the q-dependence of the matrix elements is a desirable feature 
in connection with the problem of correctly including short range correlations in a fully relativistic formalism as the 
one developed here. The use of PS vertices, which produces baryonic matrix elements only depending on the external 
variables, is not coherent with a simple but significant model of the physical mechanism behind SRC, based on the 
simultaneous exchange of a pion plus one or more heavy mesons. The correlated Feynman amplitudes involve box (or 
more complex) diagrams and one expects that the interaction matrix elements explicitly depend on the momentum 
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involved in the corresponding loop integrals. This in turn represents a strong motivation to consider PV couplings 
as the most appropriate ones for a fully relativistic approach to the NMWD of A-hypernuclei, since they prove able 
to mimic such a physical effect. 



D. Initial- and final-state wave functions 



The main theoretical ingredients entering the relativistic amplitudes of Eqs. (7) and (15) are the vertices oper- 
ators and the initial and final baryon wave functions. Since we adopt a covariant description for the strong and 
weak interaction operators, the involved wave functions are required to be 4-spinors. Their explicit expressions are 
obtained within the framework of Dirac phenomenology in presence of scalar and vector relativistic potentials. In 
the calculations presented in this work the bound nucleon states are taken as self-consistent Dirac-Hartree solutions 
derived within a relativistic mean field approach, employing a relativistic Lagrangian containing a, u, and p mesons 
contributions [68-72]. Slight modifications also permit to adapt such an approach to the determination of the initial 
A wave function and binding energy. The explicit form of the bound-state wave functions reads 



where the 2-components spin-orbital is written as 



g nK (\r\) y£{Sl) 
if nK (\r\)y^(Sl) 



yj?{il) ee £ (z w 1/2 ^ | j n) Yf'Wx?, 



/2 



(16) 



(17) 



with 



and 



1 



if K > 
if K < 



(18) 



n is the radial quantum number and k determines both the total and the orbital angular momentum quantum numbers. 
The normalization of the radial wave functions is given by 



J d\r\ |r| 2 (|/„ K (|r|)| 2 + | ff „ K (|r|)| 2 ) =1. 



(19) 



The outgoing nucleons wave functions are calculated by means of the relativistic energy-dependent complex optical 
potentials of Ref. [73], which fits proton elastic-scattering data on several nuclei in an energy range up to 1040 McV. 
In the explicit construction of the ejectile states, the direct Pauli reduction method is followed. It is well known that 
a Dirac 4-spinor, commonly represented in terms of its two Pauli 2-spinor components 



i> km (r)=( t^-fl 
■" isV ; 1 Xk,m 3 (r) 



(20) 



can be written in terms of its positive energy component <j>{r) as 



(f>k,m B { r ) 

(cr-k) 



M N + E + S(\r\)-V(\r\) 



<i>k, ms {r) 



(21) 



where S(|r|) and F(|r|) are the scalar and vector potentials for the final nucleon with energy E. The upper component 
4>{r) can be related to a Schrodinger-like wave function <j>(r) by the Darwin factor £>(|r|), i.e., 



with 



D(|r|) = 



V) == vfW 4>{r) , 

MAr + .E + S(|r|)-y(|r|) 
M N + E 



(22) 



(23) 
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The two-component wave function <f>(r) is solution of a Schrodinger equation containing equivalent central and spin- 
orbit potentials, which are functions of the energy-dependent relativistic scalar and vector potentials S and V. Its 
general form is given by 



««(|r|)^ ; | fc |(fir)_ 
Imijn 

x(!m, l/2m s |j>) Yp*{n k ). (24) 



<ftfc,m» = \/ 2^ 4tti' 



E. Decay rates 

In the complete calculation of the total and partial decay rates, as well as of polarization observables, the dynamical 
information on the elementary AiV — > NN process, given by the amplitudes in Eq. (7) or Eq. (15), are included in a 
many-body calculation for nuclear structure. The weak non-mesonic total decay rate is defined as [46, 65] 



M 2 



L N 



E\E,2 2Jh + 1 

M., H {R} 
m tl m t2 

The energy-conserving delta function connects the sum of the asymptotic energies E\^ of the two outgoing nucleons, 
coming from the underlying A7V — > NN microscopic process, with the difference between the initial hypernuclcus 
mass Mh and the total energy Er of the residual (A — 2)-particle system after the decay. A sum over Er is also 
usually understood. Integration over the phase spaces of the two final nucleons is needed, since the decay rate is 
a fully inclusive observable. Moreover, the sums in Eq. (25) encode an average over the initial hypernucleus spin 
projections Mj H , where Jh is the hypernucleus total spin, a sum over all the spin and isospin quantum numbers of 
the residual (A — 2)-system, {R} = {Jr, Mr, Tr, Mt r }, as well as a sum over the spin and isospin projections of the 
two outgoing nucleons m Sl2 and m tl 2 , respectively. If we choose a reference frame in which, for instance, the £-axis 
is aligned along the momentum fei, and exploiting the energy-conservation in the delta function, the six-dimensional 
integral in Eq. 25 can be reduced to a two-dimensional integral, one over the energy of one of the two final nucleons 
and the other one over the relative angle between the momenta of the two nucleons (due to azimuthal symmetry), 
which can be performed numerically. 

The expression for the NMWD rate T nm can be decomposed into a sum over n- and p-induced decay processes 
without any interference effects, i.e., 

Tnm = y ' ^nm[ rn t N ] ~ ^n/m + ^nm ! (26) 

m tN 

where T nm [mt N ] is defined as in Eq. (25) and |.M/i| 2 is evaluated with a fixed value of the initial-nucleon 
isospin projection, m tN = 1/2 for p- induced and m tjv = —1/2 for n- induced channels. Actually, in each term 
of JZ m m \Mfi\ 2 the m tl 2 quantum numbers are fixed, so that |.M/i| 2 would involve products of the kind 

Tf\*^[mt 11 mt 21 m' tN \ T/f ^[m tll m t21 mt N \, where, in principle, also interference effects, m' tN ^ mt N , are allowed. How- 
ever, the non diagonal products with m' tN ^ m tN are necessarily zero, since if one of the two amplitudes is non-zero 
the other one must vanish as a consequence of the charge-conservation isospin factor X of Eq. (9) (same final state 
but different initial states, pA or nA). Therefore, only the diagonal terms, m' tN = m tN , contribute and without 
interferences the coherent sum over mt N becomes an incoherent one. 

The nuclear transition amplitude, from the initial hypernuclear state to the final state of an (A — 2) residual nucleus 
and the two outgoing nucleons, is defined as 

M f i = (f\M A N^NN\i) (27) 

and can be represented in terms of the elementary two-body AA^ — y NN relativistic Feynman amplitude, 7/^ -773 
which contains all the relevant information about the weak-strong dynamics driving the global decay process. The 
final A-particle state | / ) must be further specified and decomposed into products of antisymmetric two-nucleon and 
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residual (A — 2)-nucleon wave functions. An explicit decomposition for the initial hypernuclear wave function | i ) 
can be developed following the approach introduced in Ref. [46], which is based on a weak-coupling scheme, i.e., the 
isoscalar A hyperon is assumed to be in the lsi/2 ground state and it only couples to the ground-state wave function 
of the (A — l)-nucleon core. As discussed in Ref. [46], this weak-coupling approximation has been able to yield quite 
good results in hypernuclear shell- model calculations [66]. 
The final expression for M fi is 



\M fi \ 2 [m tN ] = (T R M TR l/2m tN \T H M TH 

x A ( J C T H {\ J R T R ,j N m tN 

N 

J2 ( J R M R i N mN Jc Mc ) 

m N 

x (j c M c 1/2 m A | J h M Jh ) 2 \T f A i)n \ 



(28) 



where m tN = +1/2 for T\Z and -1/2 for T^X, {J H , M Jh ,T h , M Th } are the spin-isospin quantum numbers for 
the initial hypernucleus, {In^n} are the initial-nucleon total spin and its third component, {J C ,M C } are the same 
quantum numbers for the (^4 — l)-nucleon core, and, finally, m\ is the initial A total spin projection. In Eq. (28) 
(J c Th {\ JrTr ,j'jv) are the real coefficients of fractional parentage (c.f.p.), which allow the decomposition of the 
initial ( A — l)-nucleon core wave functions in terms of states involving a single nucleon coupled to a residual (^4 — 2)- 
nucleon state. The factor A is produced by the combination of initial- and final-state antisymmetrization factors 
with the number of AN pairs contributing to the total decay rate. Eq. (28) neglects possible quantum interference 
effects between different values of jn (and m^v), namely we are ruling out interferences between different shells (si/2 
and P3/2) for the initial nucleon. Thus the calculation does not require the c.f.p., but only the spectroscopic factors 
S = A (c.f.p) 2 , that can be taken, e.g. from Ref. [46]. 



F. Antisymmetrization and isospin factors 

A crucial role in determining the T n /T p ratio is played by the isospin content of the model, namely the X factors 
defined in Eq. (9) in terms of the SU(2) isospin operators (generally represented by the 2x2 Pauli matrices) and of 
the corresponding isospin 2-spinors for the initial A and N as well as for the two final nucleons. 

Taking advantage of the AI =1/2 isospin selection rule, from the isospin point of view, the A behaves like a neutron 
state. We can then explicitly represent the isospin spinors for the p, n and A baryons as 

x P =x7;r 1/2 = (l), xn-xr / r _1/2 = (;) ) xA= Xn . w 

With these definitions, the X isospin factors can be evaluated for all the possible combinations of the m tjv , m tl , and 
m t2 isospin projection quantum numbers. They are non-zero only for those processes in which charge is conserved, 
namely Ap — > np and An — > nn. We obtain 

X [m tN = 1/2, m tl = -1/2, m t2 = 1/2] = X^ np = -1 , (30) 

X [m tN = 1/2, m tl = 1/2, m t2 = -1/2] = X^ np = 2 , (31) 

I [m tN = -1/2, m tl = -1/2, m t . 2 = -1/2] = X^ nn = X^ nn = 1 . (32) 

All the others possibilities imply charge violation and give zero. The (d) and (e) apices refer to the direct and exchange 
diagrams of the relativistic Feynman amplitudes for the elementary processes. 

The use of the isospin formalism means that we are treating the neutron and the proton as two indistinguishable 
particles; therefore the final state is composed of two identical particles and this requires the antisymmetrization of 
the Tfl v amplitude. The antisymmetrization acts on the two final nucleons, exchanging their spin-isospin quantum 
numbers, m Si , m ti , and their momenta ki within the matrix elements defining the 7/i j7 r complex amplitude. We can 
thus define 
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where Tff\ is the Feynman amplitude for the direct diagram, given by Eqs. (5) or (12), while Tf?\ represents the 
Feynman amplitude for the exchange diagram, obtained from the same Eqs. (5) or (12), but with the interchanges 
m Sl <-» m S2 , m tl «-» m t2 and k\ o fe 2 . In addition, the antisymmetrization involves different X factors for the direct 



and exchange diagrams: = — ^ anc ^ 

Taking advantage of the factorization T\ 
as 

Ap — » np 



1 for a final nn pair. 



-(d,e) 



A P \ np = 2 for a final np pair, and l^l 1 

j(d,e) 7^'^, the antisymmetrized Feynman amplitudes can be written 



7T 4 



"Ap— mp ' /i, 7r 



r; 



j(e) ^-(e) 
Ap— ynp fi,ir 



(-) 



ri d) +2ri e) 



(34) 



An — > nn : 



r(e) 
L An- 



■f(e) 
*nn fi,7T 



V 



id) 



fi, k 



T 



(e) 



(35) 



We stress that many complex amplitudes with different quantum numbers contribute to the calculation of the nuclear 
transition amplitude. It is therefore difficult to make simple estimates of the final result. 



G. Asymmetries in polarized hypernuclei decay 



The main effect that is obtained with polarized hypernuclei is given by the angular asymmetry in the distribution 
of the emitted protons with respect to the direction of the hypernuclear polarization. It can be shown [46] that the 
non-mesonic partial decay rate for the proton-induced Ap — > np process can be written as 



p(p) 



2J H + 



- J2 (t(Jh,M Jh )=Iq(J h ), 



(36) 



where <j(Jh, Mj h ) = K / I M | i ', Jh, Mj h )\ i s the intensity of protons emitted along the quantization axis z 
for a spin projection Mj H of the hypernuclear total spin Jh- In terms of the isotropic intensity for an unpolarized 
hypernucleus, Iq(Jh), the intensity of protons emitted in the non-mesonic decay of a polarized hypernucleus (through 
the Ap —> np elementary process) along a direction forming an angle with the polarization axis is defined by 



I(9,Jh) = Wh) 1 + P„(Jh) A y (J H ) cos9 



(37) 



where P v (Jh) is the hypernuclear polarization and A v (Jh) the hypernuclear asymmetry parameter, both depending 
on the specific hypernucleus under consideration. The asymmetry A v {Jh) is a property of the non-mesonic decay 
and it only depends on the dynamical mechanism driving the weak decay. In contrast, P v {Jh) also depends on the 
kincmatical and dynamical features of the associated production reaction. The explicit expression for A v (Jh) reads 



Ay(J H ) 



3 Em Jh °{Jh,M Jh ) M Jh 



Jh + 1 Em, v(Jh,M Jh ) 



(38) 



in terms of the quantities <j(Jh,M Jh ) defined in Eq. (36). 

Within the framework of the shell-model weak coupling scheme, supposing that the A hyperon sits in the Is orbital 
and interacts (weakly) only with the nuclear core ground-state, angular momentum algebra can be employed to relate 
the polarization p\ of the A spin inside the hypernucleus to the hypernuclear polarization P y 



Pa{Jh) = ' 



Py{Jn) , 



if J 



H 



Jr 



1 

2 ' 



(39) 



where J c denotes the total spin of the (A — l)-nucleon core. It turns out useful to introduce an intrinsic A asymmetry 
parameter, a a, which should be independent of the considered hypernucleus, such that 



P y (J H ) A V (J H ) = pa{Jh)cla- 



(40) 



12 



The a a parameter removes the dependence on the hypernuclear spin Jh and is thus given by 

- Jh + 1 Ay(J H ) , if Jh = J c -\, 
• = { H \ (41) 

A y (J H ), if Jh = J c +2> 

and aA = for J H = 0. Therefore, in this weak-coupling picture, which is known to provide a good approximation 
for describing the ground state of A-hypernuclei and is particularly reliable for the non-mesonic decay (where nuclear 
structure details are not so important), a a can be interpreted as the intrinsic A asymmetry parameter for the elemen- 
tary reaction Ap — > np, involving the polarized hyperon inside the hypernucleus, and should no more depend on the 
particular hypernucleus considered. 



H. if -exchange 



The exchange of the pseudo-scalar K iso-doublets is believed to represent the dominant dynamical mechanism 
beyond the simple OPE, since the K meson is the lightest one (mx — 495 MeV) after the pion. It is thus useful to 
investigate the role of if-exchange in our relativistic model. Actually, such a contribution was included in all the most 
recent nonrelativistic calculations, and in turn proved useful to improve the theory-experiment matching concerning 
the T n /Tp decay ratio. The OKE process is driven by the two following Hamiltonians [65]: 



H™ K = l G F ml [*« x!"/ 1 / 2 (°k + C% 75) (^) f *£> 
for the weak (strangeness-changing) NNK vertex, and 



H 



(s) 

ANK 



= i 9 ANK *V 75 4>K * 



,(<>) 
A i 



for the strong (strangeness-conserving) ANK one. In the previous expressions, <px is the K meson field and 



X-l/2 ~ 



(42) 



(43) 



(44) 



is introduced as usual to enforce the AI =1/2 isospin rule. We use the Nijmegen value gANK = 12.0 for the strong 
ANK coupling [67], and the weak parity- violating and parity-conserving coupling constants are C 7 ^ = 0.76, C v ^ = 



-18.9, D% 



2.09, D p £ — 6.63 [17]. Using the PS form for the strong and weak (parity-conserving) vertices and 



proceeding as in the case of OPE, we obtain the analog of Eq. (7) for the OKE relativistic Feynman amplitude 

i G F ml qank (2tt) 5 {E 1 + E 2 - E A - E N ) 
1 I d 3 x f d 3 y 



'fi,K 



{b) (x) 



x A K (\x-y\)f(\x-y\) 

+ icjd 3 xj d 3 y [^ m J*)75V£L>) 
x A K (\x-y\)f{\x-y\) 



D P k 



UN 



(y) 



(45) 
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where the overall short range correlation function f(\x — y\) is given in Eq. (13) and the K propagator has the same 
structure as in Eq. (8), with a monopolar form factor J-k{q 2 ) for the baryon-baryon-if vertices (the same for the 
weak and strong vertices), with Ak = 1.2 GeV [65]. The isospin factor X is the same that already enters the OPE 
amplitude and it is given in Eq. (9), while the isospin factor IC is defined as 



/C = 



1/2 



X. 



-1/2 
1/2 



fy mt2 V Y m /' 



As before the K isospin factors are non-zero only for those processes in which charge is conserved, namely Ap 
and An — > nn. The quantities K. are: 



£ [m tN 
IC [m tN 



1/2, m tl 
1/2, m tl 



-1/2, m t2 
1/2, m t2 = 



1/2] 
-1/2] 



K (d) 
Ap-mp 



IC 



(e) 

Ap—t-np 



IC [m tN = -1/2, m tl = -1/2, m t2 = -1/2] 



— ' v An->nn 



/c 



(e) 

An->nn 



(46) 
-> np 

(47) 
(48) 
(49) 



The final-state antisymmetrization requires to consider an antisymmetrized Feynman amplitude, defined as the dif- 
ference of the direct and exchange diagram contributions, where the space-spin part of the exchange term is obtained 
from Eqs. (36) and (38) by interchanging the momenta fei^ and spin quantum numbers m Sl 2 of the outgoing nucleons 
scattering wave functions. The complete expressions can be found in Ref. (54). 

When we consider derivative PV couplings for the strong and weak parity-conserving interactions, the OKE 
Feynman amplitude is given by (see Eq. (15)) 



r; 



(PV') 



fi, K 



■ G F ml g ANK (2tt) S (E 1 + E 2 - E A - E N ) 



d 3 x / d 3 y 



d 3 g 

(2^)3 

x ^{q 2 )f{\x-y\) 



-iq{x-y) 



is 



(s) 
fei.r 



(x) 



(y) 



2M 



75 i V. 



AM N 



75 



+ 



iq(x-y) 



) 3 q 2 — + ie 



x Fk(q 2 )f{\x-y\) 



(s) 



fci ,m si 



is) 



(v) 



2M 



+ D 



K 



+ 



K 



+ 



2M N 



75^US, MN (y)' 



(50) 



where again q° must be fixed at q° = E A — E\ = E 2 — En- As in the OPE case, the use of PV couplings, with 
derivative terms acting on the K propagator, gives an expression for the amplitude in which the q-integral cannot 
be simply factorized, as instead happens when adopting PS vertices. The isospin factors are non-zero only for those 
processes in which charge is conserved, i.e., Ap — > np and An — > nn. When considering the (ir + K) combined effect, 
the non-mesonic decay rate becomes proportional to the squared modulus of the (antisymmetrized) sum of the OPE 
and OKE relativistic Feynman amplitudes, i.e., 



cx 



r; 



(PS,PV') 



T, 



(PS,PV')~ 



(51) 



thus the two contributions add coherently and ir-K interference effects may arise. 
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III. RESULTS AND DISCUSSION 

In this Section the main results of our model for the A 2 C non-mesonic decay are presented. We compare our 
relativistic finite-nucleus calculation with the results of the KEK experiments performed during the last five years. 

Numerical results obtained with different choices for some important theoretical ingredients arc compared to point 
out and investigate the role and relevance of each ingredient to the final results. Both PS and PV prescriptions 
for the vertex structure are analyzed, even if there arc theoretical reasons to prefer the modified PV' prescription, 
i.e., with derivatives acting on the propagator of the exchanged meson, as already discussed in Sec. II. The decay 
dynamics must be considered with great care. As a first step, we have included the OPE diagram, that is the simplest 
contribution to the process and is supposed to represent the bulk of the weak decay, and the OKE diagram, that is the 
most natural refinement of the model. The role of SRC is also discussed. The standard choice for many calculations 
is given by Eq. (11), that is a multiplicative function that represents a parametrization of a realistic AN correlation 
function obtained in the framework of a many-body calculation. We are aware that this choice can be suitable only 
for a nonrelativistic calculation and that its use in a relativistic calculation is not justified by rigorous theoretical 
considerations, but we adopt it as a useful starting point. A detailed analysis of the role of SRC is anyway beyond 
the scope of this paper. An analogous SRC function, sharing all the basic features of the AN initial one (locality, 
energy-independence, multiplicative action) can also be introduced to model short-range correlations between the two 
emitted nucleons (see Eq. (12)). 

A. Integrated observables 

In this Section we present our results for the total non-mesonic decay width T nm = T n + T p , the r n /T p ratio, and 
the <2a intrinsic asymmetry parameter. The numerical results for the OPE and tt + K models are displayed in Table 
I and II, respectively. As anticipated, we explore three possible choices for the short range baryon-baryon correlation 
functions: we can include both initial- and final-state SRC, i.e. we set /X]v( r ) as m Eq. (11) and f^(r) as in Eq. 
(12); or we take into account only initial AN correlations, i.e., we set fXj!f(r) again as in Eq. (11) and /jvjv( r ) — 1; or 
we completely switch off SRC, i.e., we set / A ™ (r) — f^(r) = 1. We notice that the configuration in which only final 
SRC are active, while the initial ones are switched off, is quite unnatural and not of particular interest: therefore, it 
will not be considered here. 

In Table III the most recent experimental results for T nm , T n / T Pl and a a are given to provide a reference scheme. 
We note that in our approach we neglect the two-nucleon induced decay channel, that can give a contribution of 
about 20-25% to T nm [41, 42, 88]. 

• 7r-exchange: As a first step, we have considered the OPE contribution, that is the starting point of all OME 
models, and we have evaluated the amplitude for both PS and PV couplings with the three possible SRC 
choices. The results are shown in Table I. The main difference between the PS and PV cases is that the use of 
PS vertices yields total decay rates T nm considerably higher than those obtained with PV' vertices: in the PS 
case T nm is more than two times larger than in the PV' case and, moreover, it overestimates the experimental 
result. The T n / T p ratio is less affected by the vertex choice: the ratio calculated with PV' couplings, and 
including initial and final SRC, is larger than the corresponding one with PS couplings by about 25%, but 
remains within the experimental range. The asymmetry is negative in both cases, when SRC are accounted for: 
negative and small with PS couplings (—0.080 or —0.064, depending on our choice to include or not final SRC 
besides initial ones) and larger in absolute value with PV couplings (—0.126 and —0.108, respectively). The 
effects of final SRC in a calculations where initial SRC are included, is generally small, independently of the 
vertices choice: the total decay rate is increased by no more than 5% and even smaller effects are obtained on 
the neutron-to-proton ratio. More significant changes (20-30%) are found for the asymmetry that, being close 
to 0, is more sensitive to the various contributions. If we compare the results with and without SRC, we see that 
without SRC, in the PS sector, T nm increases by about 25%, r„/r p increases by about 20%, and a\, which 
with SRC is small and negative, changes its sign, becoming positive and still remaining small in size. Similar 
results are obtained when the role of SRC is analyzed with PV' couplings: neglecting SRC, T nm and T n / T p are 
increased, with respect to the results obtained in presence of SRC, by about 15% and 10%, respectively, and a\ 
increases and approaches 0, still remaining negative. 

• (n + A')-exchange: The inclusion of the OKE contribution has significant, in some cases even sizeable effects 
on the integrated results shown in Table II. The decay rate r nm is somewhat reduced, more significantly for 
PV couplings (about 10%) than for PS ones. The PS choice overestimates the experimental total non-mesonic 
decay rate, while the PV choice is in fair agreement with the experiments, only slightly underestimating (when 
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TABLE I: Model results for F nm = T n + T p , V n /Y p , and «a for the \ 2 C hypernucleus when only the (OPE) (tt) diagram is 
considered. 



Model configuration 


r / p/ ree 

A nm / - 1 - A 






Pseudo-Scalar (PS) couplings 








7T 


2.426 


0.342 


— 0.080 


tt, no final SRC 


2.351 


0.344 


— 0.064 


tt, no SRC 


2.950 


0.413 


0.052 


Pseudo- Vector (PV) couplings 








7T 


0.995 


0.436 


-0.126 


7r, no final SRC 


0.965 


0.430 


-0.108 


tt, no SRC 


1.119 


0.469 


-0.029 



including SRC effects) the most recent results on T nm . We note that some underestimation is anyway expected 
in the present model, where the two-nucleon induced decay is neglected. The OKE mechanism produces a 
much larger effect on the T n /T p ratio, whose value is strongly enhanced, by a factor of about 2, in the PS 
case and reduced, by about 40%, in the PV' case. The difference between the values of T n /T p calculated with 
PS and PV couplings is therefore strongly enhanced. With PS couplings the ratio becomes now greater than 
0.6, independently of our choice for SRC, in strong disagreement with the most recent experimental extractions 
based on coincidence spectra. We have already discussed, however, why we believe that the PS choice is not 
very reliable. Adopting PV' couplings the reduction produced by the OKE mechanism is helpful to achieve a 
better agreement with the most recent experimental extraction of the T n /T p ratio, even though already the 
OPE results turned out to be satisfactory. We must notice, however, that in our calculation we get closer to 
the experimental range from above, while in most nonrelativistic calculations, where the OPE results arc quite 
small and the inclusion of the OKE mechanism increases the ratio, the same values are approached from below. 
The a\ parameter is somewhat reduced (in absolute value) with the PS choice, while with the PV choice 
it is significantly enhanced by the OKE mechanism, also turning from negative to positive in absence of SRC, 
slightly worsening the agreement with the experimental result (although both values are within the experimental 
range) . Also in this case the inclusion of OKE increases the difference between the results with the PS and PV 
couplings. All these considerations are generally unchanged by SRC, which are accounted for in the present 
model by a simple phenomenological correlation function. In general, the effects of SRC are similar to those 
obtained in the OPE model and this is independent of the vertices. The inclusion of final NN SRC, in addition 
to initial AN ones, only slightly increases the total decay rate (by just a few percent) while practically leaving 
the r„/r p ratio unchanged; a bigger effect can be seen on the asymmetry, in particular for PV vertices, even 
though the very small values of this observable makes it more difficult to draw conclusions on the possible role 
of the various theoretical ingredients in determining its final size. When SRC are completely neglected, in the 
PV case r nm and r„/r p are enhanced and a\ increases, up to 0.134. In the PS case T nm and r„/r p are 
again enhanced and a a approaches 0, still keeping the negative sign. 

Our OPE results arc different from the results produced by nonrelativistic models. Usual nonrelativistic calculations 
with OPE give small values of the r„/Pp ratio, in the range 0.05 -j- 0.2, and large negative values of the asymmetry 
parameter a\. Our rclativistic OPE calculation with PV' couplings, and including initial and final SRC, gives 
Tn/Tp = 0.436 and a\ = —0.126. The addition of OKE pulls the value of T n /T p down to 0.299 and a\ becomes 
—0.038. Both OPE and OPE+OKE results are in satisfactory agreement with the most recent determinations of the 
ratio, based on the analyses of coincidence spectra in KEK experiments (see Table III). 

Our calculation gives large values for the ratio also when using PS vertices (T n / T p = 0.342 with OPE and 0.678 
with tt + if -exchange, in presence of initial and final SRC), in contrast with the results of the relativistic calculation 
with PS couplings of Ref. [46] (r„/r p = 0.14 with OPE and 0.25 with tt + AT-cxchangc) . Different results are also 
obtained concerning the role of (initial) SRC, whose effect in the calculation of Ref. [46] is to give a strong reduction 
of the r„ m decay width (about a factor of 4) , while in our model a much more moderate reduction is found (and the 
further inclusion of final-state SRC does not change this picture). Our model is under many aspects similar to the 
relativistic one of [46], and initial SRC are described in the two calculations by the same correlation function (Eq. 11). 
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TABLE II: Model results for T nm = T n + F p , T n /T p , and a A for the A 2 C hypernucleus when both OPE and OKE (tt + K) 
diagrams are considered. 



Model configuration 


r / p/ ree 




CIA 


Pseudo-Scalar (PS) couplings 








7T ± K 


2.386 


0.678 


— 0.060 


iv + K, no final SRC 


2.277 


0.677 


— 0.054 


ty + K, no SRC 


2.847 


0.781 


0.002 


Pseudo- Vector (PV) couplings 








■k + K 


0.888 


0.299 


-0.038 


■K + K,no final SRC 


0.863 


0.299 


-0.014 


ty + K, no SRC 


1.080 


0.330 


0.134 



TABLE III: Recent experimental results for Y nm — F n ± T p , r n /F p and aA for the A 2 C hypernucleus. 



Exp. 


r / p/ ree 

A nm j J- A 


r n /r p ciA 


KEK 2004 (E307)[55] 


0.828 ± 0.087 


0.87 ±0.23 






(0.60 ± 0.23) 


KEK 2003 (E369) [25] 




0.51 ±0.15 


KEK 2004 (E508) [26, 48, 89] 


0.953 ± 0.032 


0.5^0.6 -0.16t°;28 


(single-nucleon spectra) 






KEK 2004 (E508) [29, 30] 




0.51 ±0.14 


(coincidence spectra) 






KEK 2004 (E508) [29, 30] 




0.38 ±0.14 


(coincidence spectra, [33]) 




(0.29 ±0.14) 


KEK 2004 (E508) [26] 




0.88 ±0.16 


(single-nucleon spectra, [38]) 




(0.95 ±0.21) 


Exp KEK 2004 (E508) [29, 30] 




0.46 ±0.09 


(coincidence spectra, [38]) 




(0.43 ±0.10) 



The differences between the results of the two relativistic calculations are not clear and deserve further investigation. 

The general message that can be extracted from various nonrelativistic calculations is that a good agreement with 
experiments can be reached only by considering the full OME potential and many other effects, like final state Nucleon- 
Nucleon interactions, rescattering, and intranuclear cascade, while OPE results are always far from the experimental 
values. In contrast, our calculations seem to point out a different result and suggest the relative relevance of the OPE 
and OKE contributions to get closer to the experimental results. It is not easy to explain the origin of this effect 
within our fully relativistic model. Actually, unlike what happens with the well known nonrelativistic potentials, 
where the contributions from various space-spin-isospin decay channels can be clearly isolated and separately studied, 
with our model, where we adopt a completely different approach based on the calculation of relativistic Feynman 
diagrams, a similar analysis is not possible, or, at least, this cannot be done in a simple way. 

It is not easy to establish a criterion to compare relativistic and nonrelativistic models and explain their different 
results. It would be desirable to perform a careful comparison between our relativistic description of the current 
densities and the nonrelativistic operators. We could start from our 4x4 Dirac operators and reduce them to 
effective 2x2 nonrelativistic operators. However, we believe that this would give a better understanding of our 
relativistic calculation and of its global internal coherence, but not a nonrelativistic reduction to be compared with 
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the corresponding results of nonrelativistic finite-nuclei calculations available in the literature. 

Besides final NN short-range correlations (mainly acting in proximity of the elementary interaction vertex), in 
our model only the effects of FSI due to the interaction of each one of the two outgoing nucleons with the residual 
nucleus are included, through a relativistic energy-dependent complex optical potential fitted to elastic proton-nucleus 
scattering data. The use of the same phenomenological optical potentials to describe FSI has been very successful 
in describing exclusive (e,e'p), (j,p) data and neutrino-nucleus scattering [74-87]. The imaginary part of the optical 
potential produces an absorption of flux, that is correct for an exclusive reaction, but it would be incorrect for an 
inclusive reaction, where all the channels contribute and the total flux must be conserved. An approach where FSI are 
treated in inclusive reactions by means of a complex optical potential and the total flux is conserved is presented in 
[84, 85] . In the present model the hypernuclear non-mesonic weak decay is treated as a semi- inclusive process where 
the two outogoing nucleons are detected and most of the reaction channels that are responsible for the imaginary 
part of the optical potential do not contribute. Rescattering effects where the outgoing nucleons interact with other 
nucleons in their way out of the nucleus and generate secondary nucleons can also affect the decay width. These 
processes, which are not included in the present calculations, are accounted for in [38, 49] by the intranuclear cascade 
model of [53, 54]. 

In the simplest approach, FSI are neglected and the plane- wave limit is considered for the scattering wave functions. 
In the plane-wave limit, calculations with PV' vertices give T n + T p = 2.428, T n /T p = 0.556, and aj\ = —0.210 with 
only OPE, and T n +T p = 2.263, r„/r p = 0.375, and = —0.031 with n + K exchange. It is clear that, as it 
was expected, T nm is much higher than in presence of a complex absorptive optical potential. The r„/r p ratio also 
increases with respect to the results obtained including FSI effects. This is mainly due to the fact that, though basically 
isospin-independent, these optical potentials also include Coulomb correction terms that distinguish between protons 
and neutrons. Due to the relatively small energies of the process, these terms are comparable to the central and spin- 
orbit Schroedinger-equivalent potentials and they can thus play an important role within the FSI implementation. 
The asymmetry parameter is somewhat reduced in both cases in the plane-wave limit. 

B. Kinetic energy spectra 

Further insight into the theoretical content of the model is provided by the investigation of the calculated kinetic 
energy spectra. Some examples are presented and discussed in this Section. In Figs. 2-3 the double-differential 
spectra, dT p , n /dE d9, for the Ap — >• np and An — > nn decay channels at different angles, arc displayed as a function of 
the proton and neutron kinetic energy T Pi „ = E. Here 9 is the relative angle between the momenta of the two outgoing 
nucleons. The energy spectra are calculated with PV' vertices, including both OPE and OKE contributions as well 
as initial and final SRC. It is evident from the results shown in the figures that the energy spectra for angles 6 < 120° 
display a quite flat behaviour over the whole possible energy range (from to <~ 150 MeV) and have comparable sizes, 
while the curves obtained for 9 > 120°, i.e., nearly back-to-back angles, are clearly peaked and increase rapidly in 
size with 9. Furthermore, these high-angle distributions suggest the presence of an underlying double-peak structure, 
though not very pronounced, with a first small peak at E ~ 30-^40 MeV and a much larger peak at about E ~ 70-^90 
MeV. Such a double-peak behaviour can be understood in terms of the contributions coming from the S1/2 and P3/2 
shells for the initial proton, which tend to produce different energy distributions. In the case of proton emission, the 
spectra are non-symmetric with respect to E = T max /2, where T max is the available kinetic energy for the emitted 
nucleons, due to the distorting energy-dependent optical potential acting on the final nucleons, which distinguishes 
between proton and neutron states. On the contrary, in the case of neutron-induced decay the symmetry of the 
problem, which involves the emission of two indistinguishable neutrons, leads to almost symmetric spectra. 

The kinetic energy spectra dT Pin /dE, for the proton and neutron-induced decay channels, integrated over the 
relative angle 9, are presented in Fig. 4. The results obtained with only OPE and with tt + K exchange are compared 
in the figure. The PV choice for the vertices has been adopted in the calculations, and both initial and final SRC 
are included. In the OPE case the proton spectrum is significantly larger in size than the neutron one. Moreover, the 
two-peak structure is evident in the proton spectrum, while the shape of the neutron spectrum is smoother and much 
more symmetric with respect to T max /2. The inclusion of the OKE contribution produces a more peaked proton 
spectrum but does not modify its global size. In contrast, the neutron spectrum is considerably reduced by OKE and 
its shape is slightly flattened. At the level of integrated observables, such a behaviour produces the reduction of the 
r n /r p ratio, when including kaon-exchange. In Fig. 5 the kinetic energy spectra integrated over 9 are given for PS 
couplings. The proton spectrum with only OPE again shows an asymmetric shape, with a slight two-peak behaviour 
(its global scale is now about twice the one for the corresponding PV case). The inclusion of the OKE contribution 
yields an almost identical shape but reduces the global size of the proton spectrum. The neutron energy spectrum for 
the OPE case is smoother and smaller in size with respect to the proton one. When we include the OKE contribution, 
the neutron spectrum becomes much more peaked and its size is increased. These results are in opposite trend with 
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20 40 60 80 100 120 140 

E [MeV] 



FIG. 2: Double-differential spectra (normalized to 1^ ee ) for the Ap —¥ np decay channel as a function of the kinetic energy of 
the outgoing proton, T p = E, for different values of the relative angle 9 between the momenta of the two outgoing nucleons. 
Line convention: 9 = 30° (orange), 9 = 60° (green), 9 = 90° (cyan), 9 = 120° (blue), 9 = 150° (magenta), 9 = 160° (purple), 
9 — 170° (red). Calculations are performed with PV' couplings and include n + K contributions as well as initial and final 
SRC. 




20 40 60 80 100 120 140 

E |MeV] 

FIG. 3: The same as in Fig. 2, but for neutron kinetic energy spectra, T n = E, of the An — > nn decay channel. 
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respect to what found for PV' couplings in Fig. 4 and, by consequence, in the PS case the r n /T p is considerably 
increased by the addition of the OKE mechanism and becomes larger than 0.6 (see Table II). 

r f dE MeV 

0.01 



0.008 



0.006 



0.004 



0.002 




20 40 60 80 100 120 140 

E [MeV] 



FIG. 4: Kinetic energy spectra (normalized to F^ ree ) integrated over the angle 8. Calculations are performed with PV' 
couplings and include both initial and final SRC. Red (blue) lines refer to proton (neutron) spectra, as obtained in proton- 
(neutron-)induced one-body decay channels (E = T PjIl ). Dashed lines correspond to results obtained with only OPE, solid lines 
are obtained including also OKE contributions. 

Similar analyses can also be repeated for the corresponding model configurations in which either final SRC or both 
initial and final SRC contributions are neglected. The shapes of the spectra are quite insensitive to the presence or 
absence of short range correlations, as implemented in our model. On the other hand, the global size of the spectra 
and the corresponding integrated quantities r Pi „ arc influenced, sometimes significantly, by SRC effects. 

C. Angular spectra 

In this Section we analyze the angular distributions predicted by our model. In Figs. 6-7 the double-differential 
spectra, dT p n /dE d9, for the hp — > np and An — > nn decay channels at different values of the proton and neutron 
kinetic energy, T Pi „ = E, are displayed as a function of the relative angle 9 between the momenta of the two outgoing 
nucleons. The spectra are calculated with PV' vertices, including both OPE and OKE contributions as well as initial 
and final SRC. The angular spectra for proton- and neutron-induced channels exhibit a similar behavior. In both 
figures all the curves, corresponding to different values of the kinetic energy of the outgoing nucleon, are strongly 
peaked at high angles, especially for 150° < 9 < 180°. This agrees with the fact that the elementary process driving 
the decay is a two-body A./V — > nN interaction, which preferentially yields back-to-back final nucleons. Distortion 
effects produced by the nucleon-nucleus optical potential, however, tend to smear the angular distributions, thus 
increasing the probability of emitting nucleons in non back-to-back configurations and at small angles. In the angular 
region between 0° and 120° all the spectra have similar shapes and comparable sizes, while in the back-to-back region 
they suddenly increase and differentiate among each others. The angular spectra associated with central energies (52, 
68, 84 MeV) display a higher peak, while the curves obtained for energies close to the energy endpoints (20 and 128 
MeV) are much less peaked and definitely smaller in size. 

The evidence that the most back-to-back peaked angular spectra are those pertaining to central energies can be 
understood if we observe that a proton or a neutron kinetic energy close to the middle of the available energy range 
means that the two final nucleons have approximately the same energy. Thus, due to the energy-momentum global 
conservation in the two-body process, the two outgoing particles are preferentially emitted along opposite directions. 
On the other hand, in those cases in which one of the two nucleons carries away a large part of the available energy and 
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FIG. 5: The same in Fig. 4, but for PS couplings. 
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FIG. 6: Double-differential spectra (normalized to r^ ree ) for the Ap — > np decay channel as a function of the relative angle 8 
between the momenta of the two outgoing nucleons for different values of the kinetic energy of the outgoing proton T p = E. Line 
convention: E = 20 MeV (orange), E = 36 MeV (green), E = 54 MeV (cyan), E = 68 MeV (blue), E = 84 MeV (magenta), 
E = 100 MeV (purple), E = 116 MeV (red), E = 128 MeV (yellow). Calculations are performed with PV' couplings and 
include n + K contributions as well as initial and final SRC. 
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FIG. 7: The same as in Fig. 6, but for neutron angular spectra of the An — > nn decay channel. 



the other one carries the remaining small amount, the angular correlation is weaker and the corresponding angular 
spectra are natter. The angular distributions are smeared by FSI effects. The complex energy dependent optical 
potential distorts the wave functions of the outgoing nucleons and an imbalance between the energies of the two final 
nucleons also favours a weakening of their angular correlation and a correspondingly stronger smearing of the relative 
angle distribution. 

The angular spectra integrated over the kinetic energy of the outgoing proton (for the Ap — > np channel) and 
neutron (for the An — > nn channel), dr Pi „/d0, are shown in Fig. 8. The results obtained with OPE and it + K 
exchange are compared in the figure. Calculations have been performed adopting the PV choice for the vertices 
and include both initial and final SRC. By considering only the OPE mechanism we see that the proton and neutron 
angular distributions are quite similar in shape, especially in the back-to-back region, namely from 9 = 120° — 130° 
till 9 = 180°, where they present similar peaks. In the region where 9 < 120° the angular spectra are instead 
approximately flat and strongly reduced in size. The global size of the proton spectrum is about twice the neutron 
spectrum one, coherently with the obtained results for the ratio of the integrated decay rates, T n / T p . The inclusion of 
the OKE contribution acts in opposite ways on the proton and neutron spectra. The curves are practically unchanged 
in the region 9 < 120°, apart from a slight reduction of the proton spectrum. By contrast, in the back-to-back region, 
the neutron spectrum is considerably lowered and flattened while the proton spectrum is correspondingly increased, 
becoming much more peaked towards higher relative angles. This behaviour produces the reduction of the r„/r p 
ratio shown in Table II, when the OKE contribution is included. 

The angular spectra integrated over the kinetic energy of the outgoing proton and neutron calculated with PS 
couplings (and including initial and final SRC) are shown in Fig. 9. Also in this case the OPE spectra have similar 
shapes, with the neutron curve significantly smaller than the proton one in the back-to-back region. On the contrary, 
the inclusion of the OKE mechanism induces opposite effects on the angular spectra, when considering PS and PV 
vertices: in the PS case the proton distribution is reduced and the neutron one is strongly increased. Thus, the 
r„/ T p ratio is highly enhanced and becomes greater than 0.6 (see Table II). 

A similar analysis can be repeated for the model configurations in which SRC are completely neglected. The shapes 
and the relative sizes of the distributions are strictly analogous to the corresponding results obtained when SRC are 
included. The conclusion is that, within our model, SRC do not affect in a significant way the shapes and sizes of 
kinetic energy and angular spectra. This is coherent with the use of a phenomenological SRC multiplicative, local 
and energy-independent function. 
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FIG. 8: Angular spectra (normalized to r^ ree ) integrated over the nucleon kinetic energy E = TV. Calculations are perfomed 
with PV' couplings and include initial and final SRC. Red (blue) lines refer to proton (neutron) spectra, as obtained in proton- 
neutron-) induced one-body decay channels. Dashed lines correspond to results obtained by considering only OPE, solid ones 
are obtained including also OKE contributions. 



We have presented a relativistic model for the non-mesonic weak decay of the \ 2 C hypernucleus. Over the last 
years many groups have been deeply involved in studies of hypernuclear physics and developed different nonrelativistic 
models that can satisfactorily reproduce all the experimental results. Anyway, the inclusion of many theoretical 
ingredients, like the full one-meson-exchange potential and many other ones, seems mandatory. In this paper we have 
proposed a first attempt to explain, at least in a qualitative way, all the features of the hypernuclear non-mesonic 
weak decay with a fully relativistic treatment of the weak dynamics, based on the evaluation of Feynman diagrams 
within a covariant formalism. 

We have considered the pseudo-scalar and pseudo-vector prescriptions for the weak-strong vertices involved in the 
AN — > nN elementary process. When considering the standard form of the pseudo-vector vertex, these two different 
couplings produce the same nonrelativistic limit, which is anyway inconsistent with the standard nonrelativistic 
approaches; only the modified pseudo- vector vertex, with a 4-derivative operating on the meson propagator, allows 
to obtain, in the nonrelativistic limit, the usual one-pion-exchange potential commonly employed in nonrelativistic 
calculations. In view of such considerations, we tested our model employing both pseudo-scalar and modified derivative 
pseudo-vector couplings. In our model short-range correlations are taken into account in a phenomcnological way, 
through a multiplicative local and energy-independent function. The involved wave functions are 4-spinors obtained 
within the framework of Dirac phenomenology in presence of scalar and vector potentials. Final-state interactions are 
included in the model by accounting for the interaction of each one of the two outgoing nucleons with the residual 
nucleus, that is implemented by means of a complex relativistic nucleon-nucleus optical potential. The main effect of 
the optical potentials is to produce a damping of the kinetic energy and angular spectra and to smear the two-body 
reaction kinematical correlations. 

Great care has been devoted to the decay dynamics. The pseudo-vector couplings yield predictions in reason- 
able agreement with the non-mesonic weak decay experimental rates, whereas visible discrepancies appear when the 
pseudo-scalar coupling is chosen. The role of one-pion-exchange and one-kaon-exchange diagrams has been carefully 
investigated. When using pseudo- vector couplings, the inclusion of the if-exchange is helpful in view of a comparison 
with the experimental determinations of the r„/r p ratio. On the contrary, if we adopt the pseudo-scalar vertices, the 
Zf-exchange gives puzzling big values. 

The role of initial short-range correlations is only moderate and especially visible in the total decay rates, which 
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FIG. 9: The same in Fig. 8, but for PS couplings. 



are reduced by about 20-25%, whereas r„/r p is much less sensitive to such a theoretical ingredient. The additional 
consideration of final SRC does not seem to introduce significant modifications, when initial SRC are already taken 
into account. We acknowledge that the particular implementation of SRC here adopted, namely in terms of local 
multiplicative functions inspired by non-relativistic calculations, could be unsuitable in the context of a fully relativistic 
calculation. Our choice to include such ingredients has been motivated by the relevance usually attributed to these 
correlations in non-relativistic hypernuclear decay calculations; the great care here devoted to the problem of selecting 
the right covariant vertices structure in combination with these non-relativistic SRC functions was actually aimed at 
minimizing the impact of the theoretical uncertainty introduced by such phenomenological inputs. Clearly, modeling 
short range baryon-baryon correlations directly in the relativistic framework, e.g. by considering box diagrams 
involving heavy mesons (typically the w), would definitely be a better strategy, though much more demanding from 
a calculational viewpoint: we believe this topic deserves deeper investigation and our calculation can be a starting 
point in view of similar generalizations. 

In contrast with the predictions of nonrelativistic calculations, our model produces significantly high values of the 
r„/r p ratio, especially when we only include the one-pion-exchange contribution. As a consequence, our results for 
r„/r p , both for 7r-and (jr + K )-exchange contributions, are close to the experimental measurements, without any 
apparent need to include other pseudo-scalar and vector mesons, that are usually accounted for by the full OME 
models, or to resort to more refined FSI models. 

It is not easy to understand why the obtained results are so different from well-established non-relativistic predic- 
tions. Due to the high energies involved and short distances probed, the role of relativity could really be important 
in hypernuclear non-mesonic weak decay; unfortunately, no work to date could demonstrate this point with certainty. 
A direct comparison between our model and standard nonrelativistic calculations is beyond the scope of the present 
investigation, also because it would be very difficult to establish the basis for a direct and unambiguous comparison. 
A nonrelativistic reduction would imply first of all to drop the lower components of the Dirac spinors and apply 
the proper normalization. In addition, the relativistic energy E should be put equal to M . However, this is by no 
means a nonrelativistic reduction, as relativity is directly included in the vertices and the propagators of the Feynman 
diagrams, and the nuclear current operators still involve the Dirac scalar and vector potentials, which can produce 
large differences between the results. 

We are aware that, at the present stage of development of our work, we cannot derive any definite conclusion about 
the relevance or the usefulness of a fully relativistic formalism to describe the short-range strong-weak dynamical 
mechanism driving hypernuclear non-mesonic decay. We are considering the opportunity to better describe the final 
state of the decay process, by refining the treatment of short-range strong correlations between the two outgoing 
nucleons, here simply described by means of a phenomenological non-relativistic function, and thus evaluating a 
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globally correlated relativistic wave function for the final nucleon pair. Another possible improvement is related to 
the dynamics of the model, and requires the implementation of the exchange of vector mesons, thus exploring the 
effects of these additional contributions on the integrated observables and on the simulated spectra. Anyway, we think 
that the results of our model represent an additional source of information and a partly new theoretical perspective, 
which may deserve attention and further investigation. 
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